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Abstract: We investigate the interaction of highly focused Hnearly 
polarized optical beams with a metal knife-edge both theoretically and 
experimentally. A high numerical aperture objective focus ses beams of 
various wavelengths onto samples of different sub-wavelength thicknesses 
made of several opaque and pure materials. The standard evaluation of 
the experimental data shows material and sample dependent spatial shifts 
of the reconstructed intensity distribution, where the orientation of the 
electric field with respect to the edge plays an important role. A deeper 
understanding of the interaction between the knife-edge and the incoming 
highly focused beam is gained in our theoretical model by considering 
eigenmodes of the metal-insulator-metal structure. We achieve good qualita- 
tive agreement of our numerical simulations with the experimental findings. 
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1. Introduction 

The growing interest in highly focused optical vector beams requires a proper treatment of the 
polarization state of the beam, which strongly influences the size of the focal spot |[T|. In par- 
ticular, the role of azimuthal and radial polarization has been investigated both theoretically 
||2l and experimentally (TM] and optimization strategies were proposed thereafter L5J. Further- 
more, when focusing with high numerical aperture (NA) objectives, the symmetry of the focal 
spot is broken ||3] for a linearly polarized beam and strong longitudinal components appear HJ. 
Therefore, a precise characterization of tightly focused laser beams is not just a challenge but 
is essential for further applications. In the literature many methods for beam characterization 
are described as e.g. the knife-edge method @|7], a point scan method jSl or a slit method 
ll9l [T0l etc. The data provided by the knife-edge method can be evaluated after the experiment 
in various ways as by employing an inversion algorithm involving a linear least-square method 
to measure the beam's diameter IfTTI . by performing a numerical differentiation of data in order 
to directly obtain the beam profile |[T2l or by using a direct fit with an error function llT3l . Also 
other numerical approaches lfT4l are used to characterize the beam. The knife-edge method is 
also applied to measure mechanical displacements in the nanometer range (15]. The applica- 
tion of an optical beam profiler without any moving parts using liquid-crystal displays is an- 



other example of the application of the knife-edge method lfT6l with a micro knife-edge scanner 
fabricated in a silicon-on-insulator substrate ifTTI . Beside a single knife-edge, also a periodical 
array of slits was studied recently, demonstrating the transmission anomalies of TM -polarized 

light in. 

The background of the conventional knife-edge method is the scalar diffraction theory, so 
the standard knife-edge method's evaluation scheme is polarization independent. Thus, theo- 
retically one dimensional beam scans by a knife-edge can be used in a variety of algorithms 
to retrieve two dimensional beam profiles with a Radon backward transform. However, when 
a highly focused two dimensional beam is profiled with a knife-edge, it is natural to ask our- 
selves, how well the conventional knife-edge method performs for vector-beams. The role of 
polarization in such beam measurements was studied theoretically for the knife-edge made from 
an ideal conductor 11 01 . The first precise measurements of the highly focused beams using the 
knife-edge method were performed by R. Dorn et al [3 1. However, this work already shows a 
very first experimental indication that without careful optimization of material and knife-edge 
parameters the knife-edge method can be polarization sensitive and the conventional evalua- 
tion may fail. Thus, motivated by this work we performed a systematical study on how careful 
various parameters of knife-edges for beam reconstruction have to be selected. 

The aim of our paper is to investigate in detail the interaction of highly focused linearly polar- 
ized beams with a knife-edge, made from a variety of pure materials. We investigate the beam 
profiling situations for two polarizations when the electric field is either parallel (p-polarization) 
or normal (s-polarization) to the edge. We extend previous theoretical studies on diffraction 
through a finite slit 1 19| and incorporate plasmonic modes ll20l into the model to simulate our 
experimental findings. Thus, our theory is based on the modal analysis in a metal-insulator- 
metal waveguide, which was studied in detail elsewhere 112111221 . As a result, we obtain a good 
agreement with our experiments. 

The structure of our paper is as follows. We start with the description of our experimental 
setup and explain the principle of our measurement. After that, we proceed with the discussion 
of our experimental findings. In the third chapter we discuss the physical mechanisms in the 
realizable knife-edge method and their effects on the outcome of measurements. Finally, we 
develop a theoretical model and present results of numerical simulations, which are compared 
with experimental findings. 

2. The setup, principle of measurement and experimental results 

2.1. Setup 

Most of the experiments were performed at wavelengths from 500 up to 700 nm using a tun- 
able femtosecond laser system from TOPTICA. Additionally, we performed measurements at a 
wavelength of 780 nm using a laser diode-based cw-system. The collimated linearly polarized 
Gaussian (TEMqo) laser beam was focused onto the sample using a microscope objective with 
NA of 0.9. The full width at the half maximum (FWHM) of the intensity of the incoming beam 
was 3.1 mm filling 86% of the entrance pupil of the microscope objective. The sample was 
mounted onto a piezo stage to control its 3D-position with nanometer accuracy. 

As sensors we used in-house fabricated p-i-n photodiodes. Opaque and pure material films 
were thermally deposited on a photodiode. Therefore periodic stripe-like structures were pat- 
terned using electron-beam lithography, where the slit and the stripe width is both equal (see 
Fig. [Hi- The thickness h was approx. 130 nm. For the Au samples we additionally patterned 
structures with thicknesses of approx. h ~ 70, 100, 190 nm. Forming structures directly on the 
detector surface allows the detection of a large solid angle in transmission. The materials of 
the opaque films were gold (Au), titanium (Ti) and nickel (Ni). To obtain the geometric param- 
eters of the structures after finishing all measurements we cut each structure (knife-edge) by 




Fig. 1. Electron-micrographs of one of the gold samples investigated in the experiments. 
The knife-edge width and film thickness were determined by performing cuts with a fo- 
cused ion beam (FIB) machine. The width of the investigated knife-edge is a! = 2.0 /xm, 
the slit width is / = 2.0 flm. 




Fig. 2. Schematic depiction of the knife-edge method for a two-dimensional beam (a,b). 
Typical beam profiling data (c) and their derivatives (d). The state of polarization always 
refers to the orientation of the electric field. 



using the focused-ion-beam (FIB) technique at three positions and average values of its width 
at the bottom (d = 2.005 flm) and slope angle (a = 15°) were estimated (see Fig.[T]i- The film 
thickness was measured with high accuracy by means of an atomic force microscope (AFM). 

2.2. Principle of the measurement 

The principle of the measurement is depicted in Fig. |2] The experiments are performed using 
a highly focused linearly polarized TEMoo-mode. We investigate two polarization directions 
of the incoming beam relative to the knife-edge (in the x-y plane). At first the electric field is 
oriented perpendicularly (s-polarization) and then parallel (p-polarization) to the wall of the 
knife-edge (see Fig|2](a), (b)). The investigated laser beam is blocked stepwise by the edge of 
an opaque metal stripe that is building a single knife-edge. The photocurrent generated inside 
the photodiode is recorded for each sample position xq (see Fig.|2](c)). It is proportional to the 
power P of the diffracted field detected by the photodiode 

.0 

P= / 5,(x,xo)dx, (1) 
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Fig. 3. Distance between the peaks rf,, dp versus wavelength A derived from the experimen- 
tal data for the s- and p-polarizations. The actual knife-edge width measured by SEM is 
shown by the gray bar (a). Difference of the peak positions of the photocurrent's derivative 
ds — dp versus wavelength X for various Au samples (b). Difference of the peak positions 
ds — dp of the photocurrent's derivative versus sample height h for Au samples at various 
wavelengths (c). Difference of the peak positions ds — dp of the photocurrent's derivative 
versus wavelength A for Ti and Ni samples (d). The colored lines represent results of nu- 
merical simulations. The dash-dotted lines in (b) represent a situation with noise artificially 
added to the photocurrent and smoothed with the same filter as in the experiment. The 
colored points represent experimental results. 



where is the z-component of the Poynting vector of the field at the photodiode. In the con- 
ventional knife-edge method the derivative of the photocurrent curve with respect to the sample 
position X() (see Fig.|2](d)) reconstructs a beam projection on one axis li6j|7J, so the width and 
position of the projection can be determined. 

Since the measurements with s- and p-polarized beams are performed one after another, 
a thermal drift of the sample between measurements for each polarization could introduce a 
systematic error. Therefore we profile the beam by two adjacent knife-edges, so the effects of 
sample drifts between measurements with different polarizations cancel out (see Fig. [2] (c)). 
After that the intensity of the beam projection is reconstructed from the recorded photocurrent 
and both the position and the width of the beam are evaluated. Two parameters ds and dp 
define the distance between the peaks of the reconstructed beam profiles (see Fig.|2](d)). In the 
conventional knife-edge method no polarization effects are present, so the distance between the 
peaks for s- and p-polarized beam projections is equal (ds = dp). Thus, a non zero value for 
ds — dp indicates the presence of the polarization effect. 

Measurements for each polarization are performed at 22 different positions of the knife-edge. 
Each scan line in ;ic-direction consists of 400 steps with a step size of 10 nm. For each position 
of the edge 10 data points are measured and averaged to improve the signal-to-noise ratio. The 
photocurrent data was additionally filtered using a Savitzky-Golay smoothing algorithm (11 
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Fig. 4. Difference in the peak positions — dp of the photocuiTent's derivative vs wave- 
length A for various Au samples of 130 (a) and 70 (b) nm thickness. (T) means tempered 
samples (30 sec. 400C°). (S) means a sample with a standing-alone structure to investi- 
gate a possible influence of the opposite wall on the knife-edge measurement. In this case 
other structures are missing (a,b). Theoretical difference of the peak positions of the pho- 
tocurrent's derivative (d^ — dp ) versus slit width I for Au sample (A = 130 nm) at various 
wavelengths (c). 



points) before the photocurrent curve was differentiated. 

Most scans are performed at knife-edges with a width of approx. 2 /im. Measurements on 
wider knife-edges (approximately 3 and 4 /xm) do not result in any quantitative change within 
the experimental accuracy. Hence, in what follows we solely discuss results obtained for the 2 
pirn structures. 

During the measurement the relative position of the sample is automatically calibrated every 
15 min by the determination of the position of a 500 nm hole. The hole is etched into one of 
the metal structures which is situated some microns away from the knife-edges. The detection 
of the hole position is estimated by the maximum of transmission through it. Thus long term 
stability is obtained and the actual sample position is always determined with a precision better 
than 50 nm. 

2.3. Experimental results 

We start our discussion with experimental results for a Au knife-edge with a thickness (or 
height) of h= 130 nm. The beam profiling results ds and dp are presented in Fig. |3](a) for 
two polarization states {s and p) and six particular wavelengths. In the conventional knife-edge 
theory, the zero crossing of the second derivative of the photocurrent curve coincides with the 
position X in the photocurrent curve, where half of the maximum photocurrent is reached. It 
physically indicates the position of the knife's edge EE!. So, the total distances between the 
peaks of the reconstructed profiles d^ and dp can be interpreted intuitively as the width of the 
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Fig. 5. Dependence of the ratios Wp/X (a), w^/X (b) on the wavelength A for Au samples 
of varying thickness. Dependence of the ratios Ws/X, Ws/X on the on the wavelength X for 
Ni (c) and Ti (d) samples of comparable thickness of the opaque film h= 130 nm. The col- 
ored curves represent results of numerical simulations, the colored points - of experimental 
measurements. Parameters of the model are WQp/X = 0.9h'o,5/A = 0.55. The dashed lines 
represent the FWHM's of the squared electric field estimated from the Debye integrals 1123 1 . 



single stripe if no polarization effects are present. 

As one can see in Fig. |3](a), the optical beam profiling reveals the polarization and wave- 
length dependence of the derived values ds and dp. It should be noted that the obtained knife- 
edge width for the p-polarized beam remains almost constant as the wavelength changes, al- 
though it is actually about 100 nm larger than the value determined by the SEM. Concerning 
s-polarization the knife-edge width ds depends on the wavelength and monotonically decreases 
as the wavelength increases. To characterize the polarization effect on the measurement we plot 
the value c/, — dp which also depends on the wavelength (see Fig.[3](b)). 

Experimental data were obtained for Au samples of different thicknesses. The results for the 
relative shifts ds — dp for Au are presented in Fig. |3](b) and (c). The profiling results of Au 
samples for different thickness shows a zero shift ds — dp at approx. A = 550 nm for /i = 100 
and 130 nm and also for the whole wavelength range investigated at h = 70 nm. Additionally, 
samples made from Ni and Ti were investigated. The results for Ni and Ti are very similar to 
each other (see Fig.|3](d)). As compared with Au of h = 130 nm (see Fig.[3](b)) no zero shift 
was observed for Ni and Ti samples in the investigated wavelength range. Thus, the polarization 
effect depends not only on the thickness but also on the material of the knife-edge. Due to the 
roughness of the knife-edge and the photodiode surface (see Fig.[T](b)) the method is position 
dependent, what is the main source of systematical errors. The theoretical curves, which were 
obtained using theory developed in the next part of the paper, are also presented in the Fig. |3] 
(b)-(d). Slight discrepancies between the experimental and theoretical results are due to the fact 
that in the experiment the edge position was determined as a point of maximum derivative of the 
probe current. Furthermore, the sample in the experiment had a non-zero slope angle (see Fig.[T] 



(b)). Our numerical simulations reveal a slight asymmetry in the reconstructed beam profile (see 
discussion in the next section). So, the position of the maximum of the photocurrent derivative 
does not coincide with the position of the half level of the photocurrent curve. However, if many 
measurements are averaged to reduce the influence of noise the curves become more inversion 
symmetric. Points of maximum derivative and those of half transmission start to merge. In 
fact, if we add noise to the numerical results and average afterwards the agreement between 
experimental and numerical results is even improved (see Fig|3](b)). We see a good qualitative 
agreement of our experimental and numerical findings for all samples. Only at a wavelength 
of X = 532 (see Fig.[3](c)) and a knife-edge thickness of h ~ 100 nm and 130 nm we observe 
stronger discrepancies. 

To check the reproducibility of the results several Au samples were produced. The relative 
shifts obtained in the measurements are presented in Fig.|4](a) for /i = 130 nm (sample No. 1-7) 
and (b) for li = 70 nm (sample No. 8, 9). The results on sample No. 1 were already shown in 
Fig.[3](b). Samples No. 1-3 were fabricated using the same fabrication steps. Furthermore we 
investigated the dependency of the measured relative shift on the recrystallization of the knife- 
edge material (Au) by additionally tempering (for 30 s. at 400 C°) the sample after fabrication. 
By tempering the grainy metal surface becomes smooth and the surface roughness decreases. 
The results for the tempered Au samples are shown in Fig. |4](a) (sample No. 4-6). In Fig.|4] 
(a) one can see that by tempering the knife-edges the relative shifts (t/, — dp) derived from the 
measurements are slightly decreased for short wavelengths, although, the overall qualitative 
wavelength dependency is preserved (compare samples 4-6 with 1-3). The agreement between 
the experiment and the theory is slightly better for the tempered sample of the thickness /i = 1 30 
nm, compare with Fig.[3](b). 

The aforementioned samples have a periodical geometry (see Fig.[T](a)). Therefore, we also 
check the dependence of the experimental results on the presence of an adjacent knife-edge. 
For this purpose single non-tempered knife-edges were fabricated (sample No. 7 (/j = 130 nm), 
see Fig. |4] (a) and No. 9 {li = 70 nm), see Fig|4] (b)). The measurements on these samples 
were performed under comparable conditions as described before. The results for these samples 
are similar to the results obtained with non-tempered periodically placed knife-edges and the 
differences were found to be in the range of the experimental errors. Such outcome is expected 
in our theoretical model, see next Section. The theoretical curves reveal changes of the same 
degree as in the experiment (see Figj4](c)). 

As one can see measurement results depend on many parameters including the surface qual- 
ity and roughness of the knife-edge. Therefore it is important to perform measurements of the 
relative shift at different positions of the knife-edge, which was done in all measurements per- 
formed here. 

Following our original intention we now derive the beam width from our experimental data 
(see Fig. Is). We explicitly do this for the sample No. 1 (Au, h=130 nm). To this end we mea- 
sure the FWHM of the first derivative of the photo-current signal peak (see Fig.|2](d)). As the 
beam-width is expected to depend on the wavelength linearly, we investigate the ratios Ws/X 
and Wp/X, which should be constant. As we see from the experimental data, for samples of 
thickness /i = 130 nm, the measured FWHM remains indeed nearly constant for p-polarization 
(see Fig.|5](a)) but decreases for s -polarization with growing wavelength (see Fig.|5](b)). The 
theoretical predictions are also shown in the pictures. With increasing thickness h of the knife- 
edge the width of the reconstructed profile increases (see Fig. |5]l. This accounts to the theo- 
retical prediction, that the slit acts as a complex spatial frequency filter which damps modes 
with higher spatial frequency. In order to accurately measure the beam, the spatial spectrum 
of the signal has to be preserved, while the beam is blocked by the knife-edge. Obviously this 
can be accomplished by a knife-edge of small enough thickness. Next, as the wavelength in- 



creases, the reconstructed width of the beam slowly approaches its theoretical value. However, 
the reconstruction of the p-polarized beam reveals, that the reconstructed width of the beam 
projection remains larger than its real value over the investigated wavelength range. Due to the 
rather complicated interplay of the underlying effects, we can only suggest, that this is caused 
by the remaining power-flow through the wall of the knife (see next section for more details). 
Nevertheless the agreement between the experiment and the simulations is good, however, the 
FWHM values obtained from the Debye integrals 1231 and reconstructed experimentally are 
different (see Fig. |5]l. In what follows we develop a theoretical model to gain some insight into 
the complex dynamics of light-matter interaction at a knife-edge. 

3. Theoretical model 

3.1. Realistic knife-edge versus perfect knife-edge 

We start the theoretical discussion with a brief reminder of the knife-edge method basics. In 
the original work 1;6| the following assumptions were made: the incident beam is paraxial, the 
knife-edge is made from a perfect conductor and no losses are present. For sake of simplicity we 
consider a two-dimensional situation. The incident field in the spectral domain is represented as 
plane waves with amplitudes ^(^.t) traveling at different angles a = ^xcsmkx / k. Here k = (o/cq 
is the wave vector, with k^ and k^ being the transverse and longitudinal components of the wave 
vector, CO is the frequency and cq is the speed of light in vacuum. In the first part of Figure |6](a) 
a single plane wave component k = {k^, k^) = A; (± sin a, — cos a) of the spatial spectrum S [kx] 
of the Gaussian beam is shown. Three possible interaction scenarios are shown in different 
colors. First, the plane wave part, represented by the red rays is blocked by the knife-edge 
and reflected from the metal surface. Second, the partial ray shown in green is not affected by 
the knife-edge and is properly detected. Furthermore, the wave part shown in blue experiences 
reflection (r = 1) from the wall of the knife-edge, however its amplitude remains the same and 
is therefore also detected properly. 

Next, we introduce a finite conductivity of the knife-edge but still consider paraxial beams, 
i.e. reflection r « 1 (see Fig. |6] (a)). At the side wall the boundary conditions for s- and p- 

(2) (2) 

polarizations are ~ £Ex and Ey = Ey , where superscript 2 denotes the electric field com- 
ponent in the knife-edge and e is a dielectric constant of the knife-edge material. The fields in 
the knife-edge are decaying exponentially as exp{—Kkx), where K ~ Ime'/^. They still con- 
tribute to the total detected power T as an additional term exp (^~^) ^§r^ where do is the 

FWHM of the intensity \E\^ of the scalar beam, i) = e^' for a TM incoming field and u = 1 
for the TE field. The differentiation of the transmission gives us 



and can be interpreted as the beam projection (compare with Section 2). The second term leads 
to an asymmetric profile and results in a shift of the projection maximum, which therefore 
does not longer correspond to the half maximum of the transmission curve T. The parameter 
V depends on the polarization relative to the edge. Thus, the shift of the maximum in the 
projection curve is polarization dependent in this simple model, which leads, in general, to 
a non-zero relative shift ds — dp. 

As a next step towards a realistic knife-edge, we introduce an angle dependent reflection 
r{kx) from the side wall of the knife-edge (see Fig.|6](a)). For s- and p-polarization the ampli- 
tudes of the reflected plane wave parts are Eref = Ei^crexp (i5), here 5 is an angle-dependent 
relative phase of the reflected plane wave. In the paraxial case it holds; r w 1 — Ca and 5 w Da, 




(2) 




Fig. 6. Schematic depiction of a single plane wave component k = {±kx,k-) impinging on 
a knife-edge. Here ki = {kx,k^), k2 = {—kx,kz) and r = {x,z) (a). Sketch of the considered 
structure (b). 



where C and D are different for s- (_L) and p- (||) polarization. C and D are obtained from 
Fresnel's equations and kx = ka. Thus, the reflected part of the plane wave can be written as 
Eref = exp(— ifcva + lnr — i5). Next we use Inr w — Ca. The single plane wave component 
affected by the reflection can be written as exp (ikxa) + exp [— (ikx + C + iD) a] . In the perfect 
case (Fig.|6](a)) the spatial spectra of the signal blocked by the knife-edge at the photodiode 
is Ssig (ka) = S {ka)exp{—ikxQa). However due to the reflection, the plane wave component 
—kx is replaced by the aforementioned expression. After some mathematical operations it can 
be revealed that the reflection from the wall of the knife-edge introduces further modifications 
to the signal integrated by the detector The signal experiences a polarization dependent shift 
Ax = —D/{2k), it is spatially filtered and the shape of the spectral profile is deformed. 

We note, that a possible solution to make the method less polarization dependent is the usage 
of a knife-edge with a rough surface, so the reflection will be polarization independent. For 
instance, alloys with unstable binary phases (Zn and Au) were used as knife-edge materials in 
previous investigations successfully |[3][4l. 

3.2. Field representations and eigenvalue problem 

The theoretical consideration of the previous subsection, though it is helpful for the understand- 
ing of the underlying principles, has a number of limitations. In this section we proceed with 
the development of an exact theoretical model, which considers the experimental situation ex- 
actly; vector fields and plasmonic modes, reflection from the photodiode and exact boundary 
conditions will be also included. We start with a formulation of an eigenvalue problem. 

We consider the structure represented in Fig. |6](b), where two adjacent knife-edges build a 
wide (/ > A) slit. The space is divided into three regions, of which region I (half space z > 0) 
and III (half space z < —h) are assumed to be dielectric and homogeneous with dielectric con- 
stants £i and £3, respectively. The intermediate region II (0 > z > ~h) of thickness h consists 
of an opaque (i.e. the skin-effect depth is smaller than the height h) material (for x < and 
x> I) with a dielectric constant £2 and empty space (for Z > x > 0), with the same dielectric 
constant as region I, which we put to unity £1 = 1. The considered structure is piece-wise ho- 
mogeneous and in each homogeneous part of the structure solutions of the Helmholtz equation 
describe the field propagation. A continuous solution up to first order of the derivatives is de- 
rived by combining solutions of homogeneous parts using appropriate boundary conditions. In 
the general case, in the two-dimensional waveguide structure with ohmic losses in the walls the 
coupling between transverse electric (TE) and transverse magnetic (TM) modes was reported 
II24I . So, as a further simplification, we restrict our consideration to a planar approximation, 
assuming that the incident field does not change in y direction, so the Helmholtz equation in 



each homogeneous part can be written in it's two-dimensional form 



e denotes the respective dielectric constant. The solution of equation (|3]l consists of two inde- 
pendent classes: transverse electric and transverse magnetic modes, propagating in z direction. 

We start our consideration with TE modes of the Region II. The solution of the two- 
dimensional Helmholtz equation (O for the electric field E in this case has only one nonva- 
nishing component of the electric field Ey parallel (p-polarization) to the slit walls 

BfCt) icn dEy 

H, = H- = — E, = E-= Hy = 0. (4) 

(OZq " (OZq dx 

Here Pe is the effective propagation constant of the electric field mode in the empty part of 
region II (for the sake of brevity we will call it slit) and Zq is the vacuum impedance. 

Next, the transverse magnetic (TM) modes of the metal slit have one nonvanishing compo- 
nent of the magnetic field Hy and electric field components are perpendicular (s-polarization) 
to the walls of the knife-edges 

-Hy, E, = , H^ = H-= Ey = 0. (5) 



COEoE (jc) ■ ' ' COEoE (x) dx 

Here j3,„ is the propagation constant of the magnetic field mode in the slit and e{x) = 1 for 
< X < / and e{x) = E2 elsewhere. 

In general, the eigenmodes of the metal-insulator-metal waveguide can be divided into the 
two classes: the localized and nonlocalized solutions of Eq.[3] see Ref. ||2T|| . The first class has 
a discrete spectrum of the eigenvalues and represents the field localized in the slit. The second 
class has a continuous spectrum of the eigenvalues and is necessary to ensure the boundary 
conditions on the metal surface x ^ (0,1). 

However, when light penetrates the metal knife directly from top, the electric field will decay 
exponentially in the metal ll22ll . Therefore, the non-localized solutions can be neglected when 
calculating the transmission. We verified the validity of this further simplification numerically. 
Hence, we construct only localized eigenmodes of the three-piece region II, which consists of 
three functions: the first and the last one describes plane waves propagating into the wall, the 
second one represents an interference of two plane waves inside the slit. So the eigenfunctions 
of the Helmholtz equation Ey^ and Hy y have the following form in region II 



{ Hyl (x) } ^ \ ^2^""' + ^' = 7^^^^' ^ ^ 

where v > is an integer number, enumerating the eigenfunctions. 75 = 72 + ^7? transverse 
part of the complex wave vector in the slit wall, 7i = + iyj' is the transverse wave vector in 
the slit and I is the slit width. Here, plus and minus signs account for forward and backward 
propagating waves in the slit. The constants Ci, C2, C3 and C4 in Eq. ^ are eliminated by 
enforcing the continuity of either the Ey and H^ components for the TE polarization or the 
Hy and E- components for the TM polarization. This procedure results in two characteristic 
equations 

.2r:' = ( ^l^VforTE, .^r.' = ( ^l^^VforTM (7) 



from which the modal propagation constants jSg and j3,„ are found correspondingly, ET\ . The 
electric field y and the magnetic field Hy y of the eigenfunction in TE and TM cases are 
expressed as Tiy (x)exp (±i)3(g ,„)z), where 

TTv (x) = (l - (-1) e^-\ if X e [-00,0] 
TTvW =e''i-^-(-l)V-''i(-^"-'\ ifxe [0,/] 

TTv (x) = - (-1)^) e-^^-*-^-'), if X e [Z,oo] , (8) 

where the mode index v starts with v = for TM and v = 1 for TE waves. We see that the 
mathematical expressions for the fields in TE and TM cases are functionally the same with 
the differences being hidden in the eigenvalues yi and Yi, which are derived from different 
equations, see We note here, that the boundary conditions at the walls are satisfied by Eq. 
([8]l, so Eq. (|6]l represents a field structure in Region II. The effects considered in the previous 
subsection are automatically included into the exact model. 

In general, the eigenvalues j3(g „,) of equation (|7]i are complex for complex dielectric con- 
stants. Therefore, propagating modes are attenuated in the slit. However, the eigenmodes with 
the modal number v > Vc, , where Vc- is estimated from the relation k~er — n^V^^/l'^ < 0, are 
attenuated much faster (also for lossless materials). The last note here is that solutions of eigen- 
value equation ^ in the TM case are divided into two families. The first family consists of two 
TM eigenmodes describing symmetric and antisymmetric plasmonic modes propagating along 
the walls of the slit and having evanescent field components not only in the metal but also in 
the sUt ll25l l26l . The second family of modes consists of fields, which oscillate in transverse 
direction in the slit. All eigenfunctions are orthogonal in the sense 

Ey{x).vHj^{y),^dx = dviiGvv, (9) 

where 8v^ is the Kronecker delta and Gw is an overlap integral of the electric and magnetic 
transverse fields of the eigenmode dHJ, see Ref. ||22| for more details. However the function Gw 
is in general a complex-valued one, therefore, we normalize our eigenmodes by introducing a 
real-valued norma A^v, which we define as 

Nv= J v{x) Kv (x) TT* (x) c/x, ( 1 0) 

here, i)(x) = e(x)^' for TM and d(x) ~ 1 for TE polarization. An asterisk denotes a complex 
conjugation. Those expressions are readily evaluated to 

2 _ l+e2''i'-(-l)^2en'cos7"/ Ive^i' smy'yl ve^^t'-l 

Ny ; ( ^ 1 ) Ti \ ; • (11) 

Yi Yi Yi 

here, U = £2 for TM and u = 1 for TE polarization. The normalization is, in general, not 
necessary, however, this procedure ensures a slightly better convergance in our numerics. 

In what follows we will evaluate boundary conditions between the different regions I and II 
or II and III in Fourier space. Therefore we derive the Fourier transform of the eigenfunctions 
in the slit. The Fourier image of the eigenfunction TZy (x) can be written as a sum of two images. 
The image «l>v represents the field in the walls, and 'I'v - the field in the slit 
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Fig. 7. Absolute value of the electric field |E| / |Eo| distribution on one of the / = 2 wide 
slit walls for TE (a,c) and TM (b,d) radiation at A = 532 nm (a,b) and 780 nm (c,d). The 
center of the incident beam is at j = 0. The Au knife-edge is situated at x < and the beam 
propagates in negative z-direction. The parameters of the model are WQp = 0.9woj = 0.55A. 
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5.5. Boundary value problem at the interfaces and energy flow between regions 

In this section we discuss the boundary value problem at the entrance and exit interfaces of 
the slit. The eigenfunctions of the Helmholtz equation in regions I and II are those of free 
space, namely plane waves. On the other hand the eigenfunctions (|6]l of region II have discrete 
eigenvalues. We note, that for the optimal convergence to be achieved, the continuity of the 
tangential components has to be tested differently for TE and TM polarizations 1271 . However, 
for the sake of brevity, we employ here an unified and polarization-independent approach. To 
match the fields in the three regions the problem will be investigated in the Fourier space. 

We assume, that the electric (for TE case) or magnetic (for TM case) fields in regions I and 
III can be expressed as 

F'^"' (x,z) = — L / S {k„z) exp (ik^x) dkr, (13) 



V2nJ- 

here S{kx,z) is a Fourier image of the field. The slit is illuminated from the top. Therefore in 
the region above the first interface, the spatial spectra S can be expressed as consisting of two 
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Fig. 8. a) Fraction of energy transmitted by plasmons Tp/T as a function of tlie beam dis- 
placement xq for various samples of the same thickness h = 200 nm and width 1 = 2 /xm 
at A = 532 nm (solid lines) and A = 780 (dashed lines), b-d) Dependence of the angular 
spectral width Afc^ of the transmitted signal's spectra Su-ans on the beam displacement xq. 
Thickness of the Au film is h = 200 nm, the slit width is / = 2 /im, the wavelengths are 
X = 532 nm (b), 633 nm (c), 780 nm (d). The spectral FWHM is determined by fitting 

a Gaussian function 5/,> (^.v,xo) = 5oexp 21n2 (A;.^ — ^^o)^ /A^^j to the spectral distribu- 
tion Sfrans 

{k^,xo). For comparison the black line represents the same dependence for an 
perfect knife-edge. Parameters of the theory are vvQp = 0.9woj 



:0.55A. 



parts. One part describes the field propagating into the slit 5,„, the second part the field reflected 
from the slit Syef 



S {k^,z) = Sin {kx) exp [-ij3i [k^) z] + S^ef (^.v) exp [i)3i (k^) z] 



(14) 



Here j3i = y^Sik'^ — k^ is the longitudinal component of the wave vector and Re)3 > 0, ImjS > 0. 
We assume an incoming wave with Gaussian distribution of the y component in the transverse 
plane. So at the entrance the field Ey (Hy) is 



FG(x,z = 0)=exp -{x-xoY/Wf 



(15) 



here Wo is the beam waist radius and xq is the displacement from the center of the coordinates, 
which in our case coincides with the left wall. The FWHM wo can be found as wo = V21n2Wo. 
We note, that the FWHM w of the total electric field for the TM polarization is larger than wo 
due to the presence of the component. The spatial Fourier spectra of the field (fTSt is 



Si„ (k,) = ^ exp (-i^o) exp {~klwi/4) . 
In region III light only propagates from the slit into the substrate Strc 



(16) 



S {k^,z) Straus {kx) exp [-ijSs {kx)z\ , 



(17) 



here jSs = ^Jeik^-kl. 

The electric and magnetic fields in region II consist of forward and backward propagating 
eigenmodes Ky. Thus, the resulting field is represented by the sum 



F" {x,z) = £ [flve-'^^-' + /7ve'^'' Tty (x) =Y.Fy (z) tTv (x) , (18) 



v=l 



v=l 



with the expansion coefficients fly, which describe a field propagating down to the substrate, 
whereas the by represent a reflected field. 

The energy conservation law can be written as 

Pbeam ^ Ptrcms + Pref + Pahs i (19) 

where Pheam is the power of the incoming beam. Ptrcms is the power transmitted into the region 
III through the slit. fVe/ is flie power reflected from the slit and Pai,s is the power losses in the 
slit. After some math the reflection R from the slit and transmission T through the slit can be 
expressed in terms of the spatial spectra of incoming and outgoing fields as follows 

j'_l^l5l{k,)\Sref{k,)fdk, Rc{V3)!':%j3ik,)\Srrans{kd\'dk, 

" = — I ^ 1 ^ — IT ^ J (20) 

/!»'^^ J3, (fc,) \S,„ {k^ fdk, Re(ui ) j3i (k,) \S,„ {k,)\-dk, 

where k^^ ~ eik, k^j^ = e^k. u, = 1 for TE field and u, = e^^' for TM field. The expansion 
coefficients Oy and by of the field inside the slit along with the spatial spectra Stnms, Sref are the 
unknowns. 

The problem of finding the unknowns is solved by implying the continuity of electric and 
magnetic fields at the two interfaces (z = and z = —h). Traditionally, the procedure of the 
imposition of the continuity is different for TE and TM polarizations (see Ref. 1281 l29l for 
examples). The main reasoning behind that is an optimal convergence. Nevertheless, for the 
sake of brevity, we will use the same projections for two different polarizations. This will slow 
the convergence, but a compact and polarization independent approach can be used. 

In addition to the continuity of Ey and //,,, the derivative of the electric field and the derivative 
£ ' dHy/dz of the magnetic field are continuous as well. With the help of equations (11311141) and 
dTSI l a pair of equations is obtained at the first interface (z = 0): 



1 ^ dF 

Sref (k.) = ^ (0) 0^ (^.0 + S.„ (k,) 



dF 



2^v = E i^(0)/yj+^M(0)GMv 



(21) 



where the sums run from /i = for TM and from jU = 1 for TE polarization to infinity. 

The matching of fields at the second interface with help of Eqs. iVh and ( fTSI l gives us a 
similar system of equations for the unknowns Stmns [kx) and a^, b^i 



Siram %) = exp [ift %) h] Y^Ff;L{~h)i: ^ % 



i^(-/i)7^'J-F^(-/i)G^v 



(22) 



where 

Dl U3 

4v = £° i3r ' (^-v) (^.v) Hv {k,) dk,, jfl = I" /33- 1 (^,) (^.v) Hv {k,) dk,, 

Sv = f S,„ ik,) Uy (k,) dk,, G^y = 1°° Tv {k,) Rv (^.v) (23) 

with riv defined by equation ( fT2] i. The integral Sv is an overlap integral of the incident beam 
and the eigenmode of the slit and the integral describes optical and geometrical properties 
of the slit. The overlap integral defines the behavior of the reflected and transmitted radiation. 
The explicit expressions of F^i and dFfi/dz are 

(0) = a^+b^e'^>^'\ (0) = ii3^ {-a^^b^i^^'^ 

Ff, i-h) = fl^e'^"/' + ^7^, (-/;) = ij3^ [-a,,e'l^^''+b^') . (24) 

We substitute those expressions into the second equations of the equations system (l2HI22t and 
get 

= £ e'P^^'' (g^v - j3^4 V ) + (Ciiv + jS^^li'v* ) . (25) 

Now, we rewrite Eq. (l25t in the following matrix form 

2S = N*/^a + N['*b, = N?'a + Nf ^b, (26) 
with the matrices N being defined as 

N« = diag (J3^) J' + G, n['^ = diag (e'P"'') [G - diag (j3^) J^] , 

Ni'' = diag (e'/^.M [diag (j3^) - G] , Nf ^ = -diag (j3^) - G. (27) 



They solve the resulting matrix equations for unknown vectors a and b. 

Now, from Eqs. (l2Tl i and (|22] |. we obtain the unknown Fourier images Strms and Sref- In 
particular, the transmitted power T, see (|20] |. can be expressed as 

where 

^^Mv = / (^,) n; (^,) dk,. (29) 

The physical meaning of the derived formulation can be explained as follows. The propa- 
gation constants jS^ are complex, so each eigenmode is attenuated while it propagates from 
the entrance to the exit. The incoming field excites eigenmodes in the slit. The strength of the 
excitation is defined by the overlap of the eigenmode and the incoming field given by a certain 
distribution of transverse wave vectors k^. This means, that the slit acts as a complex spatial 



filter, acting differently on the spatial frequencies according to the damping of the respective 
mode excited in the slit. The higher the value of the initial kx the stronger is the interaction of 
the light with the slit walls. Those modes are damped even stronger for deeper slits. If a certain 
thickness li is exceeded only the first mode will reach the exit of Region II, effectively leaving 
mainly two spatial frequencies kx = ii/i, kx = ±i72 in the transmitted signal (see Eq. (fT2]i). 
This holds true in general for both polarizations with one notable difference. For s-polarization 
besides the propagating modes there are two plasmonic modes, which are exponentially local- 
ized at the slit walls. For narrow slits the plasmons are damped less than all the other modes. 
Hence, the electric field distribution at the photodiode may develop two distinct spatial peaks - 
each at it's wall of the slit. 

3.4. Numerical simulations and discussion 

In the following section we present and discuss the results of our numerical simulations based 
on our theoretical model. For a given value of the beam displacement xq the procedure of 
finding the unknowns a and b (see Eq. ( |26] |) was as follows. As a first step, roots of Eqs. (|7) 
were calculated by built-in iterative MatLab procedures. The values of the propagation constant 
j3v for a perfectly conducting wall were used as an initial guess. For plasmonic modes, an 
initial guess was the propagation constant Pp = ^[£2/ + £2)] '''^ of a single plasmonic mode 
propagating at the interface between the knife-edge walls and the slit. The matrix dimensions 
were chosen in such a way that the energy conservation criterion ( fT9l ) was satisfied with a 
relative accuracy of 10^^. In general, the typical length of the unknown vectors was about 40- 
50 elements due to the beamwidth of the incident field close to the wavelength. All dielectric 
constants were taken from ll30l . We carefully choose free parameters of the theory (wop = 
0.9wo.i = 0.55A) to achieve the best possible agreement with our experimental data. The FWHM 
of the squared electric fields in the model are Wp = O.Swv and slightly differ from the expected 
in our setup ll23l . 

As a numerical result the electric field distribution at a knife-edge maid from gold is pre- 
sented in Fig.|2]for two different wavelengths and two polarizations (s and p). We remind here, 
that for small values of z the electric field distribution is approximate, see Section 3.2. The 
beam is impinging from the top and the thickness of the knife-edge is li = 200 nm. The ob- 
served patterns are similar to those observed during the diffraction of the beam from a wedge, 
see for example ||3T| . For both states of polarization beam diffraction is observed (see Fig. |2ll. 
In the case of an s-polarized incoming beam, the excitation of plasmonic modes at the edge 
is expected and can be seen in Fig. [T] (b,d). In all graphs we see a shift of the electric field 
maximum introduced by reflection from the side wall of the knife surface, as already suggested 
in the discussion in Section 3.1. The beam reaches the photodiode at different positions. In all 
cases some part of the beam penetrates the knife-edge and reaches the photodiode as well. As a 
main difference, the part of the field attributed to the plasmonic modes is largest at A = 780 nm 
for s-polarization, whereas at A = 532 nm the plasmonic modes are attenuated comparatively 
fast due to increased absorption. For that reason localization of an s-polarized field is less strong 
in the latter case. 

In addition to the Au sample we consider Ni and Ti samples of the same thickness (h = 200 
nm) and illustrate the importance of plasmonic modes in the transmission of the TM fields 
(s-polarized) at next. The amount of power transmitted by the plasmonic modes as a function 
of beam displacement is plotted in Fig. |8](a). We see that due to their spatial extent into the 
metal knife-edge plasmonic modes have a noticeable influence on the transmission, even, if 
most of the beam is blocked. So, for a gold film at A = 532 nm they account for up to 90% 
of the transmitted power for a beam impinging on the slit walls. For Ni and Ti samples this 
value decreases to nearly 70%, while at A = 780 nm the Au-knife-edge maintains it's plasmonic 



properties. The last point we have to note is that still a non-zero amount of energy is transmitted 
by plasmons even for the beam centered in the slit. This happens due to the excitation of the 
coupled plasmons pair on both walls of the slit. 

As a further result of our theoretical calculations we see that a knife edge acts as a fil- 
ter of the incoming beams spatial spectrum. To have a reference we start with a pair of 
perfect knife-edges, see Eq. for which the spatial spectrum of the transmitted signal is 
Strans {kx,xo) ~ Jq (x) exp{—ikxx)dx. For convenience, we determine the spectral FWHM 

Akx by fitting a Gaussian function Sfit {kx,xo) = S'oexp — 21n2 (k^ — k^o)^ /Ak^ to the result- 
ing spectral distribution Strans {kx,xo). The spectral profile remains symmetric for an perfect 
knife-edge and the dependence of the angular width Ak^ is presented in Fig. |8](b-d). Our in- 
tention here is to compare the spectral width of the signal chopped by the real knife-edge to 
the spectral FWHM of the same signal chopped by the perfect knife-edge. During beam pro- 
filing with a perfect knife-edge the spatial spectrum broadens, when the beam is blocked by 
the knife-edge. However, the spectral width of the signal blocked by the other knife-edges is 
slightly larger compared to the perfect case, when the beam is in the middle of the slit. Never- 
theless, as the beam approaches the knife-edge, the spectral width becomes much smaller than 
expected for a perfect knife-edge. This is a direct indication of the filtering of spatial frequen- 
cies. So, dependent on the polarization the reconstruction may result in a beam profile wider 
than expected, if other effects can be neglected. 

4. Conclusions 

In conclusion, the knife-edge method at nanoscale was investigated in detail. We experimen- 
tally as well as theoretically investigated the knife-edge technique for reconstructing the inten- 
sity profile of highly focused linearly polarized beams. Different thicknesses and edge materials 
were studied. It was demonstrated that a variety of polarization sensitive physical effects can 
influence the successful beam reconstruction: The excitation of plasmonic modes, the polariza- 
tion dependent reflection from the wall of the knife-edge and a non-zero power-flow through 
the edge material. Thus, the reconstructed beam profile strongly depends on the used material, 
the wavelength and the thickness of the knife-edge. The peak position of the conventionally 
reconstructed beam as well as its width is affected by the interaction of the knife-edge with the 
measured beam. For the fabrication of knife-edges for highly accurate beam reconstruction a 
suppression of plasmonic modes is crucial along with minimization of a polarization dependent 
reflection from the inner side of the knife-edge. For pure materials the suggested workaround 
can be the use of very thin metal films. Another possibility to minimize the influence of the plas- 
monic modes and the polarization dependent reflection is the use of films with rough internal 
structure, where the reflection is diffuse and plasmons are suppressed (compare ||3]|4l). 
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